Z + (4430) as a resonance in the Di(D[)D* channel 
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We study the possibility that the Z + (4430) state observed by Belle is a D\D* (or D[D*) resonance 
in S-wave. Focusing on its decays, We find that the open-charm decay into D*D*n is dominant. 
Furthermore, we use the re-scattering mechanism to study its hidden-charm decays and find that 
the re-scattering effects are significant in DiD* channel but not in D[D* channel. For the J P = 1~ 
candidate, with chosen parameters, we can get V(Z + — > ip' 'n + ) /T(Z + — > J/ipir + ) ~ 5.3, which tends 
to account for why the Z + (4430) is difficult to be found in J/ipn + . However, the CP candidate can 
not be ruled out by our calculations. 
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I. INTRODUCTION 

Recently, there have been a number of exciting discov- 
eries of new hadron states in the hidden-charm sector. 
Among these new states, the Z ± (4430) may be the most 
special one, which was discovered by the Belle collabo- 
ration [l[ in the invariant mass spectrum of vp'ir^ in the 
decay B -> ip'n ± K with mass m z = 4433 ± 4 ± 1 MeV 
and width T z = U±\ltu MeV. 

The distinguishable feature of Z ± (4430) is that it is 
a charged resonance-like structure. Surely it can not be 
emplaced in ordinary charmonium sector or explained by 
couple channel effects between charmonium and open- 
charm thresholds. The closeness of mz to the thresholds 
of D*(2010)Di(2420) and £>* (2010)1^ (2430) may sug- 
gest that it is induced by the threshold effects of these 
channels. However, it may not be a simple kinematic 
effect since it is very narrow and the line-shape of rp'ir 
distribution is very sharp. Presumably, it is likely to 
be a resonance of D*D\ (or D*D[). Here, a "resonance" 
means a state formed in the £>*(2010) — Di(2420) scatter- 
ing by some inter-hadron forces (say, meson exchanges), 
but not the bound state formed by the colored-confining 
forces. Note that in the B + decays, the CKM favored 
modes such as B + — >• Z>(2010)°D* + has a very large 
branching ratio of about 2%0. This can happen not only 
for the D* + but also for the excited D* + system, which 
can decay into D+K° and D°K+ as well as D+(2420)K t) 
and Di(2A2Q)K + with sizable branching ratios. Then 
in the scattering process of D*(2010)° - Z? 1 "(2420) and 
£>* (2010)° -L>? (2420), resonance-like structures Z+ and 
Z° may be formed. 

Here, we will not try to clarify the complicated dynam- 
ics in the D*D^p scattering in this paper. However, it is 

worth emphasizing that the configuration of D*D^ in 
Z should favor S-wave because a centrifugal barrier will 
prevent the formation of a resonance state. Then, the 
allowed J p of Z is CP, 1~ or 2P Thereinto, the pos- 
sibility of 2~ may be neglected since its production in 
B — > Z(4430)K decay is suppressed by the small phase 
space. 



Another unusual feature about Z(4430) is that there is 
no report [l[ on Z + — ► J/ipn + in the J/ip — > mode, 
which implies the ratio 



R 



"0' /ip 



T(Z+ -> V'7T+) 
T(Z+ -> J/V>7T+) 



» 1. 



(1) 



It is not easy to understand this " J/ -0-suppressed prob- 
lem" since J / ip has many properties similar to ip 1 and the 
phase space should favor the J/ip production. 

In this paper, we will stress on the decays of Z assum- 
ing that it is a D*D^p resonance with J p = CP or IP 

Since £)| decays to D*tt dominantly, one can expect the 
width of Z is dominated by T(Z -> D*D*n). On the 
other hand, since Z is discovered in ip'ir system with a 
quite large product branching ratio: 



B(B -> ZK)xB(Z^ ip'w) = (4.I±1.0±1.3)xlCP 5 , (2) 

the partial width T(Z — > tp'ir) should also be quite signif- 
icant. Together with the " J/ ^-suppressed problem" , this 
implies that the mechanism, which induces these hidden- 
charm decay modes, should have some nontrivial fea- 
tures. We will examine these modes in the re-scattering 
model 0, H, H, [1] , in which the D*D± component is re- 
scattered into tp^'ir by one £>(*) meson exchange. In this 
picture, the form factor suppression, which accounts for 
the off-shell effects of the intermediate meson, will favor 
ip' because its mass is closer to the D(*)D(*) thresh- 
olds than that of J/ip. We find that the lager width of 
Z — > tp'ir and the J/^-suppressed problem may be ex- 
plained if one choose adequate parameters. 

Needless to say, as charged states, the Z ± should be 
in a isospin-triplet, and their isospin partner Z° should 
exist. In our analysis, we will focus on Z + . However, all 
the results can be applied to Z° directly if the interaction 
between Z and D*D^p has isospin symmetry. 

We give the model to describe the decays of Z + (4430) 
in Sec. II. The numerical analysis and discussions are 
given in Sec. III. And a summary is given in Sec. IV. 
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II. THE DECAYS OF Z+(4430) 



As a charged resonance, the Z + (4430) can be coupled 
to either D* D^ )+ or D* + D^ )0 . We will describe their 
interactions by the following effective Lagrangians with 
isospin symmetry: 



£ ft =9 Z(£>* 



D^) + h.C. 



for CT,(3) 



e^ ff I5*I>W + 3p^ + ft-c. for 1",(4) 



where the coupling constants g 



o.i 



ZD I'D 



are blind to the 



charge of D meson, but there can be significant difference 
between these for D\ and for D[ due to the dynamics in 
D*D i p systems. As we have mentioned in Sec. I, we only 
introduce S-wave coupling between Z and D*D^ with 
J P (Z) = CP and 1~ in Eq. ^ and (g]), respectively. 
Then, the study on the decays of Z resonance is exigent 
to clarify its quantum number even its components. 

In this section, we will focus on two primary decay 
modes of Z, i.e., the open-charm mode Z — > D*D*ir and 
the hiden-charm mode Z — > J/ip(ip')ir. For the open- 
charm mode, we assume that one of the two D*s arise 
from the decay of D 1 (D[) component and the propaga- 
tion of virtual Di(D[) is described by the Breit-Wigner 
propagator. For the hiden-charm decay mode, we use re- 
scattering mechanism, which was used in the analysis of 
A(3872) -> J/iPp @. 



A. 



In the re-scattering mechanism, the decay Z + — > 
J/ip(ip')7r + can arise from exchange of a £)W meson 
between D* and Di(D[). We will not take into ac- 
count higher excited D meson medium states because 
their contributions will be suppressed by the form fac- 
tors intensively. The Feynman diagrams for Z + (4430) — > 
D* D+(D* + Dl) J/tP^')t: + are shown in Fig. [TJ The 
diagrams for D*D[ re-scattering are the same. 

To evaluate the amplitudes, we need the following ef- 
fective @, Hi Lagrangians: 



C 



-ig^D*n*{V {d^D w D*J - D*"d^D*J) 
+^D*^d^D*J ~ ^D^D*^ | ) ( 5 ) 
g.^J^-D^ +h.c, (6) 



C 



D 1 D*tv 



ijjD]>D 

ig Dl D>A-3D*Dl u d'*d' / n + (D*-D\)d 2 w 
d„D*Pd v D\ d"d v n] + h.c. 



C 



£-D*Dir 



m D *m Dl 
igD'DnD^d^irD^ + h.c. 



(8) 
(9) 



Again, we assume that the isospin symmetry is preserved 
in That is, the coupling constants are blind to the 




FIG. 1: The decay diagrams for Z+(4430) 
D*°D+(D*+D^) -> J/V#>+. 



flavor of light quarks and the coupling constants in (J7])- 
© should product by a factor of 1/V2 for interactions 
involving n°. 

For convenience, we define the following ratios: 



TD'D" 



g^'D'D* 
g^D*D* 



(10) 



In our analysis in Sec. Ill, we will not distinguish them 
seriously, but call them as "r" totally. 

All the coupling constants will be determined in the 
next section. However, it is necessary to emphasize here 
that the determinations will not account for the off-shell 
effect of the exchanged D(D*) meson, of which the virtu- 
ally can not be ignored. Such effects can be accounted 
for by introducing, e.g., the monopole Q or dipole [i[ 
form factors for off-shell vertexes. Let q denote the mo- 
mentum transferred and rrii the mass of exchanged me- 
son, the form factor can be written as 



m,i,q ) = 



A 2 -<. 



A 2 



and the cutoff A can be parameterized as 



A(mi 



oA 



QCD- 



(11) 



(12) 



As we have mentioned, this form factor may play a im- 
portant role in understanding the " J/^-suppressed prob- 
lem" . So our numerical results will be very sensitive to 
the values of n and a. Since increasing of n is equivalent 
to decreasing of a in a large kinematical region, we will 
choose to fix a — 3, and to vary the value of n from 1 to 
2. 

We are now in a position to compute the diagrams 
in Fig. [TJ If the Z+(4430) lies above the D*D l {D*D' 1 ) 
threshold, in the processes Z(px, ex) — ► D 1 (D' 1 )(pi, ei) + 
D*(p2,e2) — > ip(il}'){p3, £3) + 7r(p4), where the momenta p 
and polarization vectors e are denoted explicitly for the 
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mesons, we can calculate the absorptive part (imaginary 
part) of Fig. Q] and find it to be given by 



Abs, 



\'P2 



32ir 2 mz 
xA l (D*D 1 (D' 1 



dQAi(Z D*D 1 (D[)) 

1>W)n), (13) 



where i = (a,£>, c, d) and -pi is the 3-momentum of the 
on-shell D* meson in the rest frame of Z(4430). 

Neglecting the mass difference between charged and 
neutral D mensons, the amplitudes of Fig.QJa) and [2c) 
will be equal to those of Fig. QJb) and Hid), respec- 
tively. However, the absorptive part of the amplitude is 
indeed sensitive to the difference between the threshold 
of D*°D { ^ + and of D* + D^ )0 . This is because the ab- 
sorptive part in (|13|) is proportional to the phase space 
factor \p2\/rri2 and is sensitive to the exact position of 
the threshold. On the other hand, the absorptive part 
is strongly suppressed by this phase space factor. Even 
smeared with the Breit-Wigner distribution function 



1 



7T (t-m?) 2 + tr % (ty 



i = Di, D[, 



(14) 



the absorptive part is still less than the real part of the 
amplitude and can be neglected. 

The evaluation of the real part of the amplitude is dif- 
ficult to be achieved. We will follow Ref. [6[ to obtain it 
from absorptive part via the dispersion relation: 



Dis(m|) 



d 



a, 2 



- mi 



(15) 



where th — mi+m,2 is the threshold. To proceed, we need 
a reasonable cutoff for the integral in (fT5|) . We choose to 
add a form factor on the absorptive part 0, [l(| : 



Absj(s) — > Abs i (s)e" 



-P\?2\ 



(16) 



where the cutoff (3 can be related with the effective radius 
of interaction R by (3 = R 2 /6. The authors of Ref. [9| 
choose (3 — 0.4 GeV~ 2 , which corresponds to R « 0.3 
fm. We will choose (3 — 0.4-1.0 GeV -2 for our numerical 
analysis. 



B. Z-^D*D*tt 

For the cascade decay Z -> D*D 1 (D[) ->■ D*D*n, the 
Feynman diagrams four Z + decay modes are illustrated 
in Fig. The four diagrams in Fig. [2] are related by the 
isospin symmetry, and one can easily obtain the approx- 
imate relation: 



r(z H 



D*+D* 



T{Z+ -> D*°D* q tt + ) 



-T(Z^ 
2 V 



D*°D* + TT°). 



Since in our model, one D* directly comes from the 
virtual D^D'Y) component and have no interference with 



/ Dt / 



(a) 





(d) 



FIG. 2: The diagrams for £+(4430) -> D'D'-k^ 



the other D*, we can use the formula for cascade decay 
to evaluate the partial width as 



1 



D* + D — k+ 



(m z -m D «) 



ds^fs 



(m„+m D * ) 2 



Z+- 



( S -m 2 ) 2 + (^r w) 2 

iJ i 1 



(17) 



For the neutral state Z°, the number of the diagrams 
for decay Z° — > D*D*ir is eight. However, the total 
width of open-charm decay of Z° should be approxima- 
tion equal to that of Z + in our model. 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 

The coupling constants of DDir in©-© can be re- 
lated to the well-known parameters in Heavy Meson Chi- 
ral Lagrangian through the relations [7|, [8( : 



of / i h i 

Jit Jtt 

2 h 1 + h 2 , , . 

gD 1 D-7T = -^—n, — ^m D *m Dl , (18) 



where / w is the decay constant of tt and A x is chiral 
symmetry breaking scale. The parameters in (TT8"|) can 
be roughly estimated through the measurement of the 
width of D*+{D' 1 ,D 1 ) 2\. Using the central value of 
these widthes, the parameters \h\ + /12I, h and g can be 
evaluated and we list them in Tab. HI 

The evaluations of the coupling constants g^D* d* and 
g^ D (i) D are somewhat difficult since there are hardly any 
experimental data to be used. The two things we can use 
are symmetry and model. For g^D* d* , we can relate it to 
g-i/iD d through heavy quark symmetry, and then estimate 
its value with the help of Vector Meson Dominant (VMD) 
model [Til: 



g^D*D' ~ gipDD 



(19) 
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TABLE I: Parameters used in the calculations. 





9-ipDfD' 


\hi+h 2 \/A x 


h 


9 


8 


5 GeV 


1 Gev" 1 


0.7 


0.6 


9 ZD?D 


a ZD]'D' 


U 




r 


5 Gev 


1.5 


132 MeV 


220 MeV 


2 



where is the decay constant of ip. For g t pD 1 (D' 1 )D> 
we take it as large as 3^(44i5)DiDi which can be esti- 
mated through the prediction of r(i/>(4415) — > D\D) in 
3 P model [l|. 

The coupling constant gzD 1 D* an( i the re-scaled one 
9 ZD 1 d* I y/ mzrnr> * should be smaller than g^,u*D* since 

as a resonance in D^D*, Z must couples with D^p D* 
through some weak dynamics (say, 7r-exchange) . So we 
choose 

9 1 M = 1-5, o° f „ = 5 GeV 

ZD\' D* ' y ZD\>D* 

in Tab. U And these parameters give the prediction on 
partial widths of open-charm decays of Z(4430): 



r(z(o-) 


-> D X D* - 


■+ D*D*ir) 


= 25 MeV, 


T(Z(0-) 


-> Z^C* - 


■* D*D*n) 


= 37 MeV, 


T(Z(1-) 


-> - 


■* D*D*n) 


= 32 MeV, 


T(Z(1-) 


-> D' x £)* - 


-> D*£>*7r) 


= 46 MeV. 



(20) 



On the other hand, we also consider about that the reso- 
nance decays to D^p D* directly. After smeared with the 
Breit-Wigner distribution functions defined in (Ti"4"|) , all 
the partial widthes is about 30 GeV. So, we will use the 
numbers in (|2"D"|) in our discussions. 

As for the hidden-charm decay Z + — > ip'ir + , the width 
is squarely dependent on the value of r, which is de- 
fined in (fT0|) . It have been argued that the coupling 
between DD system and excited charmonium is gener- 
ally not weak comparing with that for J/ip. For ex- 
ample, the experimental data tell us that the coupling 
constant 5^,(377o)dd ~ 24, which is about 3 times of 
94>dd- That is, the value of r tends to be greater than 
1. In our discussion, we choose relative large value for r 
in Tab. U to examine the possibilities to understand the 
" J/"0-suppression problem" in our model. 

Evaluating the amplitudes in (|15[) numerically, we find 
that the contributions from D*D[ channel are far less 
than those from D*D\ channel. This is mainly because 
the large width of D\ reduces the threshold effects con- 
sumedly. On the other hand, the contributions arising 
from Fig. l(c)-l(d) are less than those from Fig. 1(a)- 
1(b) by a factor of 15-30. This is probably because the 
value coupling constant g^DtD which we choose in Tab. [I] 
is small. We can see it through re-scaling g^DxD by 
^Jmj^nvo and getting the number is 2.3, which is far 
less than that of g^w d* ■ 

The numerical result of the partial widths T(Z + ) — > 
i)j^-k + are listed in Tab. HT1 Here, the contributions in- 
volved are only that arising from Fig. 1(a)- 1(b) in D*Di 



TABLE II: The predictions of V{Z + ) -> Tp {,) ir + and their 
dependence on (n, jS). Here, we only involve the contributions 
arising from Fig. 1(a)- 1(b) in D* Di channel. 



(n, x GeV") 


(1.0, 1.0) 


(1.5, 0.6) 


(2.0, 0.4) 


J p (z) 


o - 


1" 


_ 


1" 


0" 


1" 


r^+(MeV) 


2.4 


4.7 


0.68 


1.4 


0.19 


0.36 


I>^+(MeV) 


1.6 


6.0 


1.1 


3.8 


0.64 


1.9 



channel. From Tab. HH one can find that the prediction 
on Rjf,/ increases with the parameter n increasing as one 
expects. For the 0~ candidate, the ratio R^/m « 3 at 
n = 2, which have implied the " J/ip suppression" in some 
degree. However, the prediction on the partial width at 
n = 2 is too small comparing with the open-charm one in 
P0|). which indicates that: B{Z + (Q~) -> - 2.5% 

and B(B -> ZK) ~ 1.6 x 10" 3 %. So the experimental 
data seems to disfavor the 0~ candidate. 

For the 1~ candidate, when n = 2, the prediction on 
the ratio R^, ^ m 5 with B(Z+((r) -► <//tt+) ~ 6%, 
which are better than that of 0~ candidate and are 
roughly consistent with experimental data. 

If we choose a < 3 [1(3] and fix n — 2, it is easy 
to get a large value prediction on R^i which can be 
consistent with the experimental constraint in (fTJ). As a 
price to pay, the prediction of T^^+ will be decreased. 
On the other hand, if the coupling constant g^DiD is not 
as small as that given in Tab. [II or if we can evaluate it 
by re-scaling g^jjMjjM y/mjj^jj, then the contributions 
arising from Fig. 1(c)- 1(d) will be comparable to those 
from Fig. l(a)-l(b), and the width r^,/ 7r + might be en- 
hanced significantly since the interference effects turn to 
be important. Moreover, in Fig. l(c)-l(d), the couplings 
of ip^DD are S-wave while those in Fig. 1(a)- 1(b) are 
P-wave. So, the configurations in Fig. 1(c)- 1(d) are more 
in favor of ^'-productions than those in Fig. 1(a)- 1(b). 

In conclusion, our calculations imply that the 
Z + (4430) state is likely to be a resonance in S-wave D*D\ 
scattering, and the spin-parity favors 1~ but can not rule 
out 0~. Further studies on its production mechanism are 
needed. 



IV. SUMMARY 

In summary, we study both the open-charm and the 
hidden-charm decays of Z(4430) based on the assump- 
tions that it is a D± D* resonance in S-wave. The par- 
tial width of open-charm decay is dominant and can be 
adjusted to be 30-40 MeV. Then, we find that the S-wave 
threshold effects arc significant in D\D* channel and are 
suppressed intensively in D[D* channel due to the larger 
width of D[. For the 1~ candidate, choosing parame- 
ters aptly, we can get T(Z+ -> ip'w+) « 2 MeV with 
R^, « 5.3, which could be roughly consistent with ex- 
perimental data. But the 0~ candidate can not be ruled 
out by our calculations. Further studies on the produc- 
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tion mechanism of Z(4430) are needed. 

Note. When this manuscript was written, two papers 
about the Z + (4430) appeared. In [l3[, a similar idea to 
ours was proposed; while in [l4j a tetraquark explanation 
was suggested. 
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